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ABSTRACT 

We apply the morphological descriptions of two-dimensional contour map, 
the so-called Minkowski functionals (the area fraction, circumference, and 
Euler characteristics), to the convergence field k(0) of the large-scale structure 
reconstructed from the shear map produced by the ray-tracing simulations. The 
perturbation theory of structure formation has suggested that the non-Gaussian 
features on the Minkowski functionals with respect to the threshold in the 
weakly nonlinear regime are induced by the three skewness parameters of k 
that are sensitive to the density parameter of matter, Q m . We show that, in 
the absence of noise due to the intrinsic ellipticities of source galaxies with 
which the perturbation theory results can be recovered, the accuracy of Q m 
determination is improved by ~ 20% using the Minkowski functionals compared 
to the conventional method of using the direct measure of skewness. 



Subject headings: cosmology: theory — cosmology: gravitational lensing 
methods: numerical 



1. Introduction 

Weak gravitational lensing caused by the large-scale structure (LSS) of the universe 
distorts the images of distant galaxies. This phenomenon is the so-called cosmic shear, 
which offers us the unique opportunity to measure directly the projected power spectrum 
of dark matter fluctuations regardless of the relation between dynamical states of the dark 
matter and luminous matter flBlandford et al. 1991| ; Miralda-Escude 1991| ; [Kaiser 1992|) . 



Recently, several independent measurements of cosmic shear have been made from deep 
'blank-field' CCD imaging surveys, and reported significant detections of shear variance 
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van Waerbeke et al. 2000| ; |Wittman et al. 2000| ; |Bacon, Rcfrcgier fc Ellis 2000fc |Kaiser 



Willson fc Luppino 2000t |Maoli et al. 2000| ). 



Due to the nonlinear evolution of density fluctuation field in the large-scale structure, 
the cosmic shear field on small angular scale is expected to display significant non- Gaussian 
features. Even in this case, for the second moment analysis it has been shown that the 
numerical results from the ray-tracing simulations are in remarkably good agreements with 
the theoretical predictions using the fitting formula for the nonlinear matter power spectrum 
( |Jain, Scljak fc White 2000| ; [Hamana, Martel fc Futamase 2000a| ). On the other hand, the 
higher order statistics can provide additional cosmological information associated with the 
non-Gaussian features. Especially, the normalized skewness parameter of the convergence 
field can be a sensitive indicator of the density parameter of matter, Q m (jBernardeau, van 



Waerbeke fc Mellier 1997| ). However, unfortunately the highly nonlinear evolution of third 



order statistics cannot be simply described by the fitting formula for the nonlinear power 
spectrum alone. Recently, the extended method that allows us to perform the skewness 
calculations in the strongly nonlinear regime has been developed using "hyper-extended 
perturbation theory" (HEPT) ([Hui 1999 ; [Scoccimarro fc Frieman 1999|) . Nevertheless, 



several works using the ray-tracing simulations have revealed that the value predicted 
by HEPT at relevant angular scales does not agree so well with the numerical results of 
skewness parameter ( |Jain, Seljak fc White 2000| ; |White fc Hu 2000| ; |Hamana et al. 2000b"D . 



Moreover, we would like to stress that it is difficult to have a physical meaning for the 
fitting formula beyond an empirical one. Therefore, it will be worth exploring again a new 
method to effectively extracting the non-Gaussian features of the convergence field in the 
weakly nonlinear regime based on the perturbation theory, which relies on a more firm 
physical basis of the structure formation. 

A possible method we propose is to use the Minkowski functionals with respect to level 
threshold; this is motivated by the fact that the functionals give the complete morphological 
descriptions of a considered field flSchmalzing fc Buchert 1997| ). For a two-dimensional 



case, the Minkowski functionals consist of the area fraction, circumference and Euler 
characteristics of the isocontour curves, where the Euler characteristics is equivalent to the 
genus statistics often used in the cosmology QGott, Melott fc Dickinson 1986|) . Recently, 



Matsubara & Jain (2000) applied the genus curve to the convergence field reconstructed 
from the ray-tracing simulations, and found that the nonlinear evolution of convergence 
induces a deviation from the specific curve of genus for the Gaussian case. On the other 
hand, the theoretical predictions based on the perturbation theory have shown that the 
non- Gaussian features on the Minkowski functionals are completely characterized by the 
skewness parameters of the convergence field in the weakly nonlinear regime ( [Matsubara 



2000| and see also equation ([!])). These results offer a possibility to extract the skewness 
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parameters using the Minkowski functionals of the reconstructed convergence field. The 
purpose of this Letter is thus to investigate how accurately Q m can be determined from the 
skewness parameters estimated by fitting the numerical results to theoretical predictions of 
the Minkowski functionals. 



2. The Ray- Tracing Simulation and the Minkowski Functionals 

We use shear and convergence fields modeled from the ray tracing simulations through 
the dark matter distribution of N-body simulations following the previous methods by 
Hamana, Martel fc Futamase 2000"a] and |White fc Hu 200Q| . The original N-body simulations 



of the large-scale structure were performed with the P M code (see |Jing fc Fang 1994 



and |Jing 199q in detail). The following discussions focus on two cosmological models, 
summarized in Table 1, and we used three different realizations for each model. As for the 
power spectrum of matter fluctuations, we assume the cold dark matter (CDM) model with 
the transfer function given by Bardeen et al. (1986) and the shape parameter T = Qoh. All 
the simulations employ 256 3 (~ 17 million) particles in a (100/i _1 Mpc) 3 comoving box and 
start at redshift z\ = 36. The gravitational softening length e is 39/i _1 kpc. 

We use the multiple lens-plane algorithm to follow the propagations of light rays 
through the simulated matter distributions. In this algorithm, the matter content of each 
box at a certain redshift is projected onto a single plane perpendicular to the line of sight. 
We use typically ~ 20 equally spaced lens-planes in the comoving distance between source 
and observer. The particle positions on each plane are interpolated onto a grid of size 
2048 2 . In order to avoid possible correlations between different lens-planes, in each plane 
we choose one of the three realizations at the considered redshift and then project the mass 
distribution along a randomly chosen one of the three coordinate axes, translate the mass 
distribution by a random vector, and randomly rotate it in a unit of 7r/2. We consider a 
set of lens-planes between the source and observer as a different realization and use ten 
such realizations to estimate the cosmic variance associated with the measurements of weak 
lensing fields. Further details of the ray-tracing simulation are given in Hamana, Martel & 
Futamase (2000). 

The fields we use are 3° on a side. Each light ray is traced by the Born approximation 
and hence can be handled as a straight line that radially extends from observer. Throughout 
this Letter, we assume that all source galaxies are at a redshift of z s = 1 and that 
their number densities is n = 30 arcmin -2 . We then make the cosmic shear field, 7(0), 
on each grid from the ray-tracing simulations, and perform the smoothing on 7(0) by 
using a top-hat filter. Using the relation between the Fourier transforms and 7(0), 
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k(l) = [(If — ^)7i(0 + Zhh^il)}/! 2 , and assuming the periodic boundary condition, k(0) 
is reconstructed on each grid from the cosmic shear field. Figure [I] shows examples of 
the reconstructed convergence field. To compute the Minkowski functionals, we label the 
convergence field by the threshold value 2/(0) that is defined by v(0) = k(0)/<jo, where <To 
is the rms of k defined by a\ = (k 2 )- 

In a two-dimensional case, the Minkowski functionals are the area fraction Vo(v), 
circumference t>i(z/), and Euler characteristics ^2(1^) for the isocontour curve with threshold 
v that fully characterize the morphology of the field. The Euler characteristic is a purely 
topological quantity, which is equal to the number of isolated high-threshold regions minus 
the number of isolated low-threshold regions. To calculate the Minkowski functionals for 
the reconstructed convergence field given as pixel data, we employed the method developed 
by Winitzki & Kosowsky (1997). 

On the other hand, under the hypothesis that the initial perturbations are Gaussian as 
supported by the inflationary scenario, Matsubara (2000) recently derived the analytical 
formula of the Minkowski functionals based on the perturbation theory that can be applied 
to the weakly nonlinear convergence field: 

v (u) = ierfc [^= \ + „ !^ e ~' y2/2 ^o g oH 2 (^) 
vi(u) = 

= ^^^^Hi^ + M^^ (!) 

where 0\ is defined by a\ = ((V/t) 2 ) and H n (z/) is the nth order Hermite 
polynomial, so, si, and s 2 denote the skewness parameters defined by so = (k 3 )/o"o, 
si = -(3/4)(k 2 (V 2 k))/(^ct 2 ), and s 2 = -3((Vk • \7 k)(V 2 k)) /a\ , respectively, where the 
quantity Sq is the skewness parameter conventionally used in the previous works of weak 
lensing. Equation (|T[) indicates that those skewness parameters can be new statistical 
indicators of the deviations from the specific Gaussian predictions of vq(v), V\{y) and 
V2(u) with so = si = S2 = 0. It should be noted that so, s± and s 2 themselves can be 
given as functions of the cosmological parameters for the CDM model and the smoothing 
scale of the top-hat filter (Pernardeau, van Waerbeke fc Mellier 1997|) , which reveals that 




the skewness parameters are particularly sensitive to Q m . Therefore, we propose that 
comparing the theoretical predictions ([!]) to their numerical (or observed) results could 



place constraints on the cosmological parameters. In some previous work (e.g., |Matsubara 



|fc Jain 2000| ), the area fraction to labeling the Minkowski functionals has been used instead 



of the threshold in order to cancel out the horizontal shift of those functionals that is due 
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to the nonlinear evolution of the underlying density fluctuations on the high threshold side. 
However, this operation merely means that the area function Vq{v) for the non-Gaussian 
field is transformed closer to its specific curve for the Gaussian case. For this reason, we do 
not employ the operation and simply use the threshold v. 

3. Results 

Figure ^| shows both the analytical and numerical results of the area fraction vq(u) 
(left panels), circumference Vi(v) (middle panels), and Euler characteristics v 2 (v) (right 
panels) per square arcmin as a function of the threshold v for the convergence fields 
with two different smoothing scales of 9 = V (upper panels) and 9 = 8' (lower panels), 
respectively. In those plots, normalizations of the analytical predictions except vq{v) are 
determined by minimizing the x 2 value for the fitting between the predictions (|l|) and the 
numerical results. The mean values and error bars in each bin of v are estimated from the 
ten different realizations with the area of 3 x 3 square degrees, and the error corresponds 
to the cosmic variance associated with the measurements of the Minkowski functionals. 
Non-Gaussian features on the functionals for the noise-free convergence field are due to 
nonlinear gravitational clustering; at negative v it has a cutoff related to the minimum 
k resulting from empty beams and has a tail at positive v due to collapsed halos. For 
the small smoothing scale of 9 = V, there are large differences between the analytical 
predictions and the numerical results. This is because the highly nonlinear evolution of 
the density field has a large effect on the convergence field. For the large smoothing scale 
of 9 = 8', on which the convergence field is expected to be in weakly nonlinear regime, 
the numerical results are broadly consistent with the analytical predictions. Note that the 
reason that the result of 9 = 8' has larger error bars than that of 9 = V is due to the fewer 
number of statistical samples. 

Figure |3| shows the values of s , S\ and s 2 calculated by the perturbation theory, 
the direct measurement of sq (top left) from the reconstructed convergence field and 
the estimations of sq (top right), s\ (bottom left) and s 2 (bottom right) obtained from 
the x 2 fitting between the theoretical predictions ([I]) of the Minkowski functionals and 
their simulation results. Here we have used only the simulation data in the range of 
— 1.5 < v < 1.5, because we expect that the convergence field in this range is still in 
the weakly nonlinear regime and therefore can be applied to the perturbation theory 
predictions. In these figures, assuming that the survey of weak lensing is performed over 
the area of 9 x 9 square degrees, we estimated the error bars by multiplying the variance 
directly obtained from the ten realizations as shown in Figure |2| by a factor of 1/3. We 
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have confirmed that the measurement of Euler characteristics, is also sensitive to the 
discreteness effect of pixel data. Therefore, to minimize the unresolved uncertainties, 
we determined the parameters of s , s±, and s 2 in the following procedure. First, we 
determine so from the fitting of Vq(v) because the non-Gaussian features of vq(v) in the 
theoretical prediction (|l|) depends on s and cro, where <To is also computed directly from 
the reconstructed convergence field according to the definition <j\ = (k 2 )- Similarly, by 
using the already determined value of Sq, we determine Si from the shape of v\. Finally, we 
use the shape of Euler characteristics t> 2 (^) to determine the s 2 parameter. Note that this 
fitting procedure causes the large error of s 2 . The top left panel in Figure |3] shows that for 
all smoothing scales the direct measurement of sq tends to largely overestimate the value of 

50 calculated by the perturbation theory. This is because the direct measurement is more 
sensitive to the strong nonlinear rare events in the convergence distribution such as halos 
of dark matter. On the other hand, for SCDM model with 9 = 2', 4' and 8', the values of 
so obtained from our method using vq(v) fairly improve the estimations for so predicted 
by the perturbation theory. For comparison, thin lines in the top left panel of Figure [| 
show the direct measurement of So in the same range of v (—1.5 < v < 1.5) as used in our 
method. It is still clear that the modified direct measurement of So also fails to predict 
its value from the perturbation theory for all the smoothing scales. Similarly, the values 
of Si obtained from our method are very similar to the values of s\ from the perturbation 
theory for the smoothing scales of 9 £ 2'. However, one can see that the result of s 2 from 
our method cannot reproduce the value of the perturbation theory mainly because of the 
fitting procedure described above, and the results of s and Si for the smallest smoothing 
scale of 9 — V do not work well. For these reasons, we will not use the results of s and 

51 for 9 = 1' and s 2 for the determination of fi m . On the other hand, it apparently seems 
that the errors of the skewness determinations for ACDM model are larger than those of 
SCDM. This result comes from the fact that the skewness variation Aso = 10 for the flat 
universe models around ACDM model corresponds to | AQ m \= 0.05, while around SCDM 
model Aso = 0.9 corresponds to the same | AQ m |. Actually, as will be shown, the relative 
accuracy of the Q m determination is not so different in both SCDM and ACDM models. 

Table ^| summarizes the results for the Q m determination with a best-fit value and la 
error, which are obtained from the direct measurements of so and from the estimations of 
So and s± using the Minkowski functionals for the smoothing scales of 9 = 2', 4' and 8', 
respectively. We here employed the current favored flat universe models with Q m + Q\ = 1. 
The table clearly shows that our method improves the accuracy of Q m determination by 
~ 20% compared to that determined from the direct measure of skewness. 
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4. Discussion 

In this Letter we addressed the issue of how accurately the density parameter, fl m , can 
be determined from the non-Gaussian signatures in the simulated weak lensing field based 
on the perturbation theory of structure formation instead of the empirical fitting formula. 
For this purpose, we have shown that the Minkowski functionals of convergence maps 
reconstructed from the cosmic shear field can be a useful new method. This is because the 
Minkowski functionals can effectively pick up the weakly nonlinear non-Gaussian features in 
the appropriate range of threshold, in which the perturbation theory can be safely applied. 
In fact, our numerical results have shown that the fl m determination of using the Minkowski 
functionals produces ~ 20% accurate best-fit value to the input value of fl m compared with 
the result of using the direct measurement of skewness. However, we still have to further 
investigate possible uncertainties due to the limited number of numerical realizations used 
in this Letter by increasing the number, and this will be our future work. 

In this Letter, we have not considered the effect of intrinsic ellipticities of source 
galaxies on our method. Nevertheless, for the practical purpose, it is critical to take into 
account this effect, and therefore we will need the theoretical predictions of the Minkowski 
functionals, including the noise effect. This study is now in progress and will be presented 
elsewhere. In practice, it will also be necessary to take into account the redshift distribution 
of source galaxies. However, previous works have quantitatively shown that, even if using 
a more realistic model for the redshift distribution of source galaxies as expressed by 
n(z) oc z 2 exp[— (z/zq) 2 ' 5 } with the mean redshift of unity, the magnitude of cosmic shear 
signal is changed only by ~ 10% compared to the result of using all the sources distributed 
at z s = 1 (e.g., |Jain, Seljak fc White 2000[) . Therefore, we prospect that the change of 
source distribution does not largely affect our results. 
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Fig. 1. — The reconstructed convergence maps used to compute the Minkowski functionals 
for a field of 3° on a side. Two cosmological models are shown: SCDM model (left) and 
ACDM model (right). 



This preprint was prepared with the AAS IATgX macros v4.0. 

Table 1: Summary of the parameters used for the N-body simulations, h is the Hubble 
constant in units of 100 kms' 1 Mpc -1 . 



Model 




Q A 


r 


o-s 


m p (/ i " 1 M )) 


SCDM 


1.0 


0.0 


0.5 


0.6 


1.7xl0 10 


ACDM 


0.3 


0.7 


0.21 


1.0 


5.0xl0 9 
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Fig. 2. — Minkowski functionals, the area fraction (left), circumference (middle), and Euler 
characteristics (right) per square arcmin for the reconstructed convergence fields with two 
different smoothing scales of 9 = V (upper) and 9 = 8' (lower), respectively. Solid (SCDM 
model) and dashed (ACDM) curves show the analytical predictions of Minkowski functionals 
calculated by equation ([I]). The mean values (circles and boxes for SCDM and ACDM 
models, respectively) and the error bars in each bin of v are estimated from the ten different 
realizations with the area of 3 x 3 square degrees. 



Table 2: The value of Vt m estimated from the direct measurement of sq and from the 
measurements of so and s\ through the fitting of Minkowski functionals for simulation data 
with 9 = 2', 4' and 8'. We here employed the flat universe models with Q m + Q\ = 1. 



Model 


fl m from the direct measure of sq 


fl m from Minkowski functionals 


SCDM (n m 


= 1.0) 


0.50 ± 0.19 


0.78 ± 0.22 


ACDM (O m 


= 0.3) 


0.24 ± 0.05 


0.31 ± 0.07 
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Fig. 3. — The directly measured value of skewness parameter (top left) from the simulated 
convergence maps, and the skewness parameters s (top right), si (bottom left), and s 2 
(bottom right) estimated from fitting between the simulated Minkowski functionals and their 
analytical predictions. In each panel, solid lines with and without circle symbols represent 
the results by simulation and by the perturbation theory, respectively, in SCDM model, 
while dashed lines with and without box symbols are results in ACDM model. The error 
bars are estimated from ten realizations for the observation field with 9° x 9° area (see text 
in detail). In addition, in top left panel thin lines represent the direct measure of so in the 
same range of —1.5 < v < 1.5 as used in the Minkowski functionals method. 



